Strain-induced topological insulator phase transition in HgSe 
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Using ah initio electronic structure calculations we investigate the change of the band structure 
and the topological invariant in HgSe (non-centrosymmetric system) under two different type 
of uniaxial strain along the [001] and [110] directions, respectively. Both compressive [001] and 
[110] strain leads to the opening of a (crystal field) band gap (with a maximum value of about 37 
meV) in the vicinity of F, and the concomitant formation of a camel-back- (inverse camel-back-) 
shape valence (conduction) band along the direction perpendicular to the strain with a minimum 
(maximum) at F. We find that the invariant vq = \ which demonstrates conclusively that HgSe 
is a strong topological insulator (TI). With further increase of the strain the band gap decreases 
vanishing at a critical strain value (which depends on the strain type) where HgSe undergoes a 
transition from a strong TI to a trivial (normal) insulator. HgSe exhibits a similar behavior under 
a tensile [110] uniaxial strain. On the other hand, HgSe remains a normal insulator by applying a 
[001] tensile uniaxial strain. Complementary electronic structure calculations of the non-polar (110) 
surface under compressive [110] tensile strain show two Dirac cones at the F point whose spin chiral 
states are associated with the top and bottom slab surfaces. 

PACS numbers: 71.15.Ap, 72.25.Dc, 73.20.At, 73.43.Nq 



I. INTRODUCTION 

There has been an explosion of theoretical and experi- 
mental research in the search for three-dimensional (3D) 
topological insulators (TIs), which are a new quantum 
state of matter with promising applications in spintronics 
and multifunctional topological-based devices and inter- 
faces pai2j. TIs are characterized by gapped bulk states, 
similar to a normal (trivial) insulators, with robust gap- 
less topologically protected metallic surface states result- 
ing from the strong spin-orbit coupling (SOC). This in 
turn leads to an odd number of Dirac-cone-like disper- 
sions crossing the Fermi level at the surface and to a 
non-trivial Z2 topological invariant. The surface states 
exhibit an unconventional spin texture with the spin di- 
rection locked to the two-dimensional (2D) electron mo- 
mentum. Current ongoing research efforts focus primar- 
ily on the prototypical single Dirac-cone family of Bi2Te3, 
Bi2Se3, and Sb2Te3 3D TIs, where angle- resolved- pho- 
toemission experiments have shown the chiral spin tex- 
ture of the surface states [3H5]- 

Although the band inversion near the time-reversal in- 
variant momenta (TRIM) points (fc = —k mod G, where 
G is a reciprocal lattice vector) , is a strong indication that 
a system is TI, it is not definitive because the topological 
invariant is a global property of the entire occupied Bril- 
louin zone (BZ). Thus, in order to confirm conclusively 
whether a 3D system is topologically trivial or non-trivial 
it is essential to compute the Z2 invariant. For band in- 
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sulators with crystal inversion symmetry, the topological 
invariant number can be easily computed as the product 
of half of the parity (Kramers pairs have identical par- 
ities) numbers for all the occupied states at the TRIM 
points[nj. For the general case of non-centrosymmetric 
systems (such as zinc-blende structures) the topological 
invariant can be determined either by integrating both 
the Berry's connection and curvature over half of the 
BZ[7 or by calculating the evolution of the Wannier func- 
tion center during a "time-reversal pumping" process [S]. 
These two equivalent approaches are discussed in detail 
in Sec. II. 

The mercury chalcogenides HgSe and HgTe in the zinc- 
blende (ZB) structure belong to a group of unique mate- 
rials exhibiting the so-called inverted band structure. In 
this band structure, the energy of the anion s-derived Fg 
state, which is the conduction-band minimum (CBM) at 
F for most ZB semiconductors, is below the anion p-like 
valence-band maximum (VBM) with Fg symmetry. How- 
ever, the experimental results for the electronic structure 
of HgSe are controversial. Early photoelectron spectro- 
scopic experiments [9, 10] on n-type HgSe indicated that 
this material should have a gap with a positive value for 
Eq = E{re) - E{rs) = 0.42 eV. On the other hand, pho- 
toemission experiments [11] of the non-polar HgSe (110) 
surface (to minimize the effects of surface charges) could 
not unequivocally determine the sign of Eq due to lack 
of resolution of the sequence of the Fg, F7 and Fg bands. 

Currently, the consensus of a large amount of exper- 
imental results [12] for HgSe seems to give the sequence 
of the energy bands in ascending order is Fg, F7 and 
Fg, with Eq « -0.2 eV and the spin-orbit splitting, 
Aso = E{Ts) - EITy) w 0.4 eV. Ab initio calculations 
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based on the GW (G is the Green's function and W is the 
screened Coulomb interaction) scheme find p"3l - fl6] Eq « 
-0.32 to -0.58 eV and Aso ~ 0.27 to 0.32 eV, in agree- 
ment with experiment. 

In contrast, in HgTe the sequence of the energy bands 
is Fy, Fg, and Fglin]. Interestingly, recent transport 
experiments f 17' on a 70-nm-thick HgTe film, strained 
by epitaxial growth on a CdTe substrate, show that the 
strain induces a band gap in the otherwise semi metallic 
HgTe, rendering it a 3D TI. Similar strain appears when 
thin HgSe films are epitaxially grown on CdTe or CdSe 
substrates. 

The objective of this work is to investigate the effects 
of uniaxial strain on the topological phase of HgSe em- 
ploying ab initio calculations. The rest of this paper is 
organized as follows: Section II describes the methodol- 
ogy used to calculate topological invariants and the band 
structures. In section III we analyze the evolution of the 
band structure and the topological invariant number un- 
der two different type of uniaxial strain along the [001] 
and [110] directions, respectively. In order to corroborate 
the topological insulator nature we also present results 
for the electronic structure of the non-polar (110) HgSe 
surface and the spin polarization of the surface states. 
Finally, conclusions are summarized in Sec. IV. 



II. METHODOLOGY 

For 2D systems, there are two commonly employed 
numerical approaches to determine whether a non- 
centrosymmetric insulator is a TI. 

1. Use the periodic part of the Bloch functions, 
\un{k)), to calculate the Berry connection, A{k) — 
i^„(u„(A:)|Vj|M„(fc)), involving the sum over the 
occupied bands. [7] For a 2D system, both the Berry 
connection and the Berry curvature, ^{k) = x 

A{k), can then be used to directly compute the 
topological invariant, z, involving the integration 
over half of the BZ, denoted by 03 i.e.. 



1 
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mod 2. 



2. Employ a recently proposed equivalent method for 
the Z2 topological invariant based on the U(2N) 
non-Abelian Berry connection[5]. This approach 
allows the identification of the topological nature of 
a general band insulator without any of the gauge- 
fixing problems that plague the concrete, previ- 
ous implementation of invariants. The main idea 
of the method is to calculate the evolution of the 
Wannier function center (WFC) directly during a 
"time-reversal pumping" process, which is a Z2 
analog to the charge polarization proposed by Fu 
and Kane [SI. The evolution of the WFC along ky 



corresponds to the phase factor, 6', of the eigen- 
values of the position operator, X, projected into 
the occupied subspace. Each state of the nth oc- 
cupied band is indexed by three quantum numbers 
|n, kx, ky) . This allows to define a square matrix 
F{kx,ky) containing the overlap integrals, 

{F{kx, ^a))mn ~ ^y\^^ ~^ Akx, ky) , 

where Akx = is the discrete spacing of Nx 
points. One can in turn calculate the complex uni- 
tary square matrix 

A'x-l 

D{ky)= n FijAkx,ky), 
3=0 

which has the complex eigenvalues \i{ky) — 
|A;(fcj,)|e'^'''^"''. The topological invariant z is then 
calculated by the even or odd number of crossings of 
any arbitrary horizontal (0=const.) reference line 
with the evolution of the Ois, mod 2. 

For a 3D bulk system, it is necessary to calculate two 
different invariants, zq and z^r, for the two different BZ 
planes, fc^ = and k^ = tt, containing eight TRIM points. 
The computation for each plane is analogous to the 2D 
case which has four TRIM points. We employ the second 
aforementioned approach to calculate the 3D topological 
invariant, which is defined as = (zq — -^tt) mod 2. 
The system is in the strong topological insulator (trivial 
insulator) phase, if vq is 1 (0). More specifically, the sys- 
tem is a non-trivial or strong TI if the evolution curves of 
the WFC cross an arbitrary reference line an odd num- 
ber of times in the A:^ = plane (zg — 1) and an even 
number of times in the kz = tt plane (z^ = 0), or vice 
versa[S]. On the other hand, if the evolution curves cross 
the arbitrary reference line an even (zo^z^r^O) or odd 
(zo=z^=l) number of times in both the k^ = and k^ = 
TT planes, the system is topologically a trivial or normal 
insulator (NI). 

The ab initio calculations employed the Vienna Ab 
initio Simulation Package (VASP) code[THl [12] employ- 
ing the Perdew-Burke-Ernzerhof (PBE) flavor of the 
exchange-correlation functional [20] and the projected- 
augmented-wave approach '2 1| to represent the electron- 
ion interaction. The spin-orbit coupling was included in 
the self-consistent calculations. The Hg 5c? electrons are 
treated explicitly as valence electrons. The energy cut- 
off of the plane-wave expansion of the basis functions 
was set to be 550 (450) eV for the bulk (surface) cal- 
culations. The BZ integration was performed with a 
11x11x11 (7x7x1) k-point Monkhorst-Pack meshp2] 
for the bulk (surface) calculations. The equilibrium vol- 
ume, Vo = 61.86A per formula unit, is in good agree- 
ment with previous ab initio calculations 23J. The ex- 
perimental lattice constant a=6.084 A 23J differs by 3% 
from the theoretical equilibrium of PBE. The non-polar 
(110) surface is modeled by a periodic slab consisting of 
21 atomic layers with a 12 A-thick vacuum region sepa- 
rating the periodic slabs. 
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FIG. 1: Band structure relative to the Fermi energy of un- 
strained HgSe along the F-X ([100]) and F-L ([111]) symme- 
try directions of the BZ [fc^ = ^^(1,1,1) and kx = ^(1,0,0), 
where a is the equilibrium lattice constant] calculated by (a) 
the PBE and (b) the MBJLDA approach [29], cmbj = 1.2. 



An accurate description of the electronic structure is a 
prerequisite in the search and discovery efforts for next- 
generation TIs. It is known that density functional the- 
ory studies [21H2B] of TIs based on local or semi-local 
exchange-correlation functionals can incorrectly predict 
a material to be a TI when in reality it is not [27] , to the 
detriment of experiment-theory interplay. This is usu- 
ally due to the fact that such a local or semi-local treat- 
ment systematically underestimates the band gap and 
can lead to an incorrect ordering of the frontier bands 
at the TRIM points, which is the determining factor in 
the prediction of the topological phase. Furthermore, 
it can also give rise to wrong band topologies and ef- 
fective masses|28j. Thus, we have also carried out cal- 
culations employing the Modified Becke- Johnson Local 
Density Approximation (MBJLDA) method ^5]. which 
predicts[28 band gaps, effective masses and, most impor- 
tantly, frontier-band ordering that are in very good agree- 
ment with the computationally more intense GW(13j and 
hybrid- functional[30j approaches. 

Figures [ija) and[l];b) display the PBE and MBJLDA 
band structures of the un-deformed HgSe, respectively, 
along two symmetry directions, where the MBJLDA pa- 
rameter c=1.2. PBE predicts the correct level ordering 
in the upper valence-band region with Eq = -1.2 eV and 
Ago = 0.25 eV, thus yielding a Fg state which is lying 
far low relative to the Fg state. On the other hand, the 
MBJLDA shifts the Fg state to a higher energy, result- 
ing in Eo = -0.23 eV and Ago — 0.20 eV, in agreement 
with the values of Eo = -0.32 eV and Aso = 0.27 eV 
of recent GW calculations[13 . The fact that both the 
PBE and the LDA 13 exchange correlation functionals 
for HgSe yield the correct frontier band ordering (Fy, 
Fg) allows for a reliable determination of the topologi- 
cal invariant [3] . Thus, throughout the remainder of this 
paper the electronic structure calculations under strain 
are carried out using the PBE exchange correlation func- 
tional. The atomic positions were fully relaxed using the 
conjugate gradient algorithm until all interatomic forces 
are smaller than 0.01 eV/nm. 



III. RESULTS AND DISCUSSION 
A. Bulk band structure under uniaxial strain 

In order to open a band gap in the semi-metallic HgSe, 
one needs to employ strain. We have carried out calcula- 
tions under hydrostatic pressure in a wide volume range, 
from -30% to -|-30% of the equilibrium volume Vq, in 
which the equilibrium structure of HgSe remains zinc- 
blende [23]. We find that within this wide range HgSe al- 
ways has a zero band gap and hence remains a semimetal. 

Therefore, we have investigated the evolution of the 
electronic structure and the topological invariant under 
uniaxial strain, e[^j./] = along the [001] and [110] 

directions which lowers the Td symmetry of the ZB struc- 
ture to D2d and G2v , respectively, and opens a band gap 
near the F point. In both cases, the in-plane lattice con- 
stant perpendicular to the strain was fixed to that of the 
equilibrium ZB structure. Our ab initio calculations of 
the total energy as a function of both uniaxial strains 
indicate that the corresponding structures are mechani- 
cally stable up to strain of ±15%. 



1. [001] strain of D2d symmetry 

Figs. [2]^a) and (b) show the band structures of HgSe 
under e[ooi] of -5% and -)-3%, respectively, perpendicular 
to the strain direction. Under a [001] uniaxial tensile 
(compressive) strain e[ooi] > (<)0, of D2d symmetry, the 
half-filled four-fold degenerate Fg state splits into an two- 
fold degenerate occupied valence Fy^ (Fg^) and a two- 
fold degenerate unoccupied conduction Fgc (Ptc) energy 
levels leading to a positive (negative) splitting, Ajp^j^j = 
E{V^)-E{T^) > (<)0. 



Under compressive strain (Ajg^j^j < 0) the maxima of 
the two highest non-degenerate valence bands along the 
[110] direction, corresponding to the irreducible represen- 
tations 3 and 4, respectively, occur away from F. This 
gives rise to a dumbbell- or camel-back-shape valence 
bands with valence band minima at F, due to the fact 
that band energies along [110] (not shown in the figure) 
are equal to those along [110] due to time reversal symme- 
try. Similarly, the S4 conduction band displays a reverse 
camel-back shape with a maximum at F and a minimum 
along the S direction away from F. These valence bands 
along S become double-degenerate at F. This leads to 
an indirect gap of Egap = 37 meV located in a direction 
perpendicular to the strain axis, in contrast to the case 
of the [001] tensile stress in Fig. |2](b), where the energy 
gap Egap occurs along the strain axis. 

Fig. ^h) shows the band structure of HgSe under ten- 
sile strain of e[ooi] — +3% (Ajg^^j > 0) along the high 
symmetry [110] and [001] directions parallel and perpen- 
dicular to the strain direction, respectively. Note that 
along the A ([001]) direction of symmetry both the 
conduction and valence bands are doubly degenerate (de- 
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FIG. 2: (Color online) Band structure of HgSe (relative to the 
Fermi energy Ep) under uniaxial strain along the [001] direc- 
tion of T)2d symmetry near the F point along the symmetry 
directions parallel and perpendicular to the strain direction, 
(a) The compressive strain of £[001] = -5% yields a TI {i^o ~ 1) 
with Egap ~ 37 meV. (b) The tensile strain of £[001] = +3% 
yields a NI {vq = 0) with Egap = 13 meV. The red and blue 
color denote the two spin polarized bands, while the green 
refers to non-spin polarized bands. 

noted by the green curves) of A5 representation. The 
crystal-field-splitting gap, Egap = 13 meV, which is al- 
lowed by the same symmetry of the bands ^Slj , is observed 
away from F. 



2. [110] strain of C2v symmetry 

Application of uniaxial stress along the [110] reduces 
the symmetry further to the orthorhombic C2V As this 
space group does not allow four-fold-generate states, the 
ZB Fg splits under strain and creates a gap between a va- 
lence T^y and conduction Fs^, both of which are two-fold 
degenerate. Under expansion the frontier energy bands 
exchange their ordering. 

Figs. [3]ja)-(d) show the evolution of the band struc- 
ture around F, under a uniaxial [110] strain of (a) —5%, 
(b) +3%, (c) -1-6% and (d) +8%, along the [110] {x + y) 

and [110] {—X + y) high k symmetry directions, which 
are parallel and perpendicular to the strain direction, re- 
spectively. Under £[110] = —5% the direct gap at F is 
Aj^j^j^Qj = E{r^c) ~ -E'(F5„) — 116 meV. The minimum 
band gap Egap = 31 meV occurs perpendicular to the 
strain direction between the camel-back-shape valence 
and inverse camel-back-shape conduction bands of S4 ir- 
reducible representation. 

Figures |4]ja) and (b) display the evolution of the band 
gap Egap and the energies of the conduction and valence 
band states i?(F„) at F as a function of uniaxial strain 
(-15% < e < +15%) applied along the [001] and [110] di- 
rections, respectively. The gray shaded areas denote the 
strong topological insulating phase with vq = I, while the 
white background corresponds to the normal insulating 
phase {vq = 0). The evolution of with uniaxial strain 
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FIG. 3: (Color online) Band structure of HgSe under uniaxial 
strain along the [110] direction of C2v symmetry near the F 
point along the high symmetry directions parallel and per- 
pendicular to the strain direction, (a) The compressive strain 
of £[110] ~ —5% yields a TI {vq = 1) with Egap = 31 meV. (b) 
The tensile strain of £[110] = +3% yields a TI {vq = 1) with 
Egap ~ 9 meV. (c) Critical tensile strain e[iio] ~ +6% where 
the conduction and valence bands of irreducible representa- 
tion 4 cross at Ef and HgSe undergoes a TI -f-^ NI transition, 
(d) At larger tensile strain of e[iio] = +8% HgSe becomes a 
NI {1^0 — 0) with a band gap Egap = 25 meV. The red and 
blue color denote the spin polarized bands. 



is discussed in more detail in the next section. For [001] 
strain, the frontier valence and conduction bands at F are 
of Fg and F7 symmetry, while they are of F5 symmetry for 
the [110] strain. Under both types of uniaxial compres- 
sion Egap exhibits a non-monotonic behavior with strain, 
reaching its maximum value in the topological non-trivial 
phase (z/q = 1) at e[ooi] = -5% and ^[iio] = -6%, and van- 
ishing at some critical strain of £[001] = -10% and e[iio] = 
-13%, respectively, where the system undergoes a transi- 
tion into the trivial insulator phase. On the other hand, 
HgSe has a different response between [001] and [110] 
tensile (e >0) strain. Namely, for e[ooi] > the band gap 
is trivial (fo = 0) and increases with strain [Fig. |4][a)] 
, while £[110] > the system is in a strong topological 
insulator phase {vq — 1) for e[iio] +5% and it becomes 
a trivial insulator {vq = 0) at larger strain [Fig. |4|^b)]. 
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FIG. 4: (Color online) Variation of the band gap Egap and 
the energies of the conduction band and valence band states, 
E{Tn) at the F point as a function of uniaxial strain along (a) 
the [001] direction of D2d symmetry (where n= 6 and 7) and 
(b) the [110] direction of C2v symmetry (where n = 5c and 
5v). The gray shaded areas denote the non-trivial topological 
phase, where vo = 1. 



B. Evolution of Topological Invariant 




FIG. 5: (Color online) Schematic of the strain-induced evo- 
lution of the band structure along an arbitrary direction 
(F ^ i^K) in the kz=Q plane where the spin-split bands are 
denoted by red and blue. (a)The strain splits the originally 
four-fold degenerate frontier bands at F in electrons and holes 
bands which does not involve any band crossing and hence the 
bulk system is in a topologically trivial phase with zq = 0. 
(b) Band inversion of one of the spin-split bands (blue) at 
the critical strain, (c) Odd number (one) of band inversions 
for one of the bands (blue) at a larger strain giving rise to 
a band gap opening away from F and the concomitant for- 
mation of a camel-back- (inverse camel-back-) shape valence 
(conduction) band along the F -o- ±K direction, rendering the 
system a strong TI (zq = 1). Note that the inverted bands 
exhibit both electron and hole character. 



An intriguing question is whether there is a correlation 
between the change of the frontier valence and conduc- 
tion bands and corresponding changes in the topological 
invariant vq. Adiabatic deformations of the band struc- 
ture (not involving band touching or overlaps across the 
band gap) under external perturbations (such as strain) 
leave the topology invariant. Murakami [32] has recently 
studied the classification of all the possible gap closing in 
2D and 3D and showed that the gap closings is associated 
with the change of the topological invariant number. 

In order to elucidate the underlying mechanism of the 
evolution of both the band inversion and i^o, we show 
schematically in Fig. [5] the strain- induced change of the 
band structure close to the TRIM T point along an ar- 
bitrary direction (T o -izK) in the fcz=0 plane, where 
the spin-split bands are denoted by red and blue. Fig. 
[sja) shows that under, for example, a uniaxial tensile 
strain, the originally four-fold degenerate Fg state (un- 
der zero strain) splits in the spin-polarized electron and 
hole bands without involving any band inversion. This 
corresponds to the band structure along the [110] direc- 
tion for £[001] > in Fig. [2|b). In order to reach the 
topological phase and its characteristic counterpropaga- 
tion of opposite spin states (i.e. — —S^) on the 

surface, it is necessary to "knot" (entangle) the frontier 
bands. [23 El] That is achieved by inverting an odd (one) 
number of spin-split bands (blue) in Fig. ^h). Further 
increase of the strain [Fig. [sjc)] gives rise to the opening 
of a (crystal field) band gap with the concomitant forma- 
tion of a camel-back- (inverse camel-back-) shape valence 
(conduction) band along the F o ±K direction. This 
corresponds to the non-degenerate E4 valence and con- 
duction bands along the [110] direction for £[110] < in 



Fig. [sj^a) . Since no other band inversions occur along the 
[110] direction, the systems is a strong topological insula- 
tor (2:0=1). Similarly, there is only one band inversion of 
the valence and conduction bands of E4 character along 
the [110] direction for ^[110] >0 in Fig. IT|b), rendering 
the system a strong TI. Moreover, the doubly-degenerate 
valence and conduction bands of A5 character along the 
[001] direction for e[ooi] > [Fig. [sFb)] have an even 
(two) number of band inversions and nence the system is 
a topologically a trivial insulator (zo=0). 

While the above mechanism, based on the changes of 
the frontier electronic bands (eigenvalues), provides an 
intuitive explanation of the topological phase transitions, 
it requires knowledge of the overlap of the 2D conduction 
and valence energy surfaces. It becomes less practical 
when these surfaces are not smooth and present multi- 
ple extrema. That is the case for Figs. [3]^b)-(d), for 
instance, requiring a careful examination along multiple 
fc-directions. For such cases, a less intuitive, but more 
stringent criterion, is the determination of the topolog- 
ical invariant, = {zq — Zt^) mod 2, based on the evo- 
lution of the Wannier functions (which include a phase 
factor), as explained in section II. Nonetheless, one notes 
a clear trend in the evolution of the bands along the [110] 
direction. The gap between the valence and conduction 
bands of irreducible representation 4 is progressively re- 
duced with increasing tensile strain until they cross at ~0 
eV for £[110] = 5 — 6% [Fig. [sj^c)]. The system undergoes 
a second band inversion that "unknots" the band struc- 
ture. As expected, this closing of the gap coincides with 
the topological phase transition predicted by the invari- 
ant vq. For larger strain, this gap increases monotonically 
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FIG. 6: (Color online) Evolution of Wannier centers for HgSe 
along ky in the kz=0 plane under uniaxial strain along the 
[001] direction of (a) e[ooi] ~ -5% and (b) epoi] = +15%, and 
along the [110] direction of (c) £[110] = -5% and (b) £[110] = 
+8%, respectively. The evolution lines cross any arbitrary 
reference line parallel to ky (for example the red dotted line) 
an odd (even) number of times yielding zq — 1 (0). In all 
cases the corresponding = 0. A^^^ — 48. 



[Fig. [Sjjd)] and without changes to the topological phase. 

In HgSe, band inversions occur only in the vicinity of F, 
because the conduction and valence states at other TRIM 
points are too far separated in energy. Consequently, the 
bands in the kz=iT plane are always not inverted [anal- 
ogous to Fig. [Sjja)] and hence -2^=0 in the entire strain 
range, thus yielding i/q = (z^^ — zq mod 2) = zq. 

Figures [6]ja) and |6][b) show the evolution lines of the 
Wannier centers along ky in the kz—Q plane under a uni- 
axial compression of ejooi] = -5% and e[ooi] — -15% along 
the [001] direction, respectively. Under the -5% (-15%) 
compression the evolution curves (blue) cross any arbi- 
trary line parallel to the horizontal axis (for example the 
red dotted line) an odd (even) number of points, thus 
yielding zq = 1 (zq = 0). Furthermore, the evolution 
lines of the Wannier centers along ky in the kz='!r plane 
under both uniaxial compression (not shown here) yield 
Zt^ = 0. This demonstrates that under uniaxial compres- 
sion along [001], there is a topological difference between 
the the TI phase (1/9=1) at e[ooi] = —5% and the trivial 
insulating phase {vq=0) at £[001] = —15%. Interestingly, 
the critical value of £[001] = —10% where HgSe under- 
goes a TI NI transition is close to the critical value 
of £ ^ —11% where the ZB structure undergoes a tran- 
sition to the cinnabar struct ure|23j. In Figures [6]jc) and 
[6f^d) we show the evolution of the Wannier centers in the 
kz=0 plane under a uniaxial compression (£[110] = —5%) 
and expansion (£[110] = +8%) along the [110] direction, 
where zq = 1 and zq = 0, respectively. Since in both 




[110]^ r ^[110] 



FIG. 7: (Color online) Band structure of the (110) HgSe 
surface under a [110] uniaxial compressive strain (£[110] = 
— 10%) along the symmetry directions in the two-dimensional 
BZ showing a single Dirac-cone at F at 52 meV below the 
Fermi energy. The gray and blue colored bands correspond 
to bulk- and surface-derived states, where the latter are asso- 
ciated with the two top and bottom atomic layers [55]. 

cases Zt^ — 0, HgSe is a TI {i/q — 1) under —5% uniaxial 
compression while it is in the trivial insulating phase {vq 
= 0) under 8% uniaxial expansion. 



C. Surface band structure and spin polarization 

While the Z2 criterion is sufficient to ascertain whether 
a 3D bulk system is a TI, an additional commonly em- 
ployed criterion is the existence of gapless surface states 
that have spin texture. Therefore, we have carried out 
electronic structure calculations of the non-polar (110) 
HgSe surface, to minimize the effects of surface charges 
and reconstruction associated with polar surfaces. The 
band structure of a uniaxially strained (£[110] — -10%) 
HgSe (110) slab is shown in Fig. [t] along the symmetry 
directions in the 2D BZ. The states marked in blue are 
identified to be spatially confined to the top and bottom 
boundary layers of the slab, i.e. surface states, using the 
approach of Park et aZ[35]. The 21 atomic layer slab ex- 
hibits a mirror symmetry with respect to both surfaces. 
One can clearly see that the topological surface states 
form two superimposed Dirac cones at the F point, each 
spatially located on the top and bottom surfaces. The 
Dirac point (52 meV below the Fermi energy) is four- 
fold degenerate due to Kramer's degeneracy at the TRIM 
point F. The inset shows the variation of the energy 
gap of the slab at F as function of the slab thickness, 
D — Na/\/2, where N is the number of atomic layers 
and a — 6.278 A is the equilibrium lattice constant. The 
large value of the energy gap in thin slabs is due to the 
interactions between the states localized at the opposite 
surfaces. As the number of layers increases, the size of 
the gap decreases converging to ~ 1 meV at the critical 
thickness of about 7 nm (A'' ~ 16). 

Fig. [8] shows constant energy contours of the surface 
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-0.1 0.1 -0.1 0.1 



FIG. 8: (Color online) Constant energy contours of surface 
spin polarization vector of the uniaxially strained (110) HgSe 
surface in Fig. [7| within an energy of ±25 meV about Ef 
(all state above the Dirac point), projected onto the (a) top 
and (b) bottom (Oil) HgSe surfaces, respectively. The red 
arrow contours denote the spin polarization vector of states 
at Ef separating those with energies higher and lower than 
Ef- The units of k are 2n/a. 

spin polarization vector projected on the — ky plane 
in the vicinity of the Fermi energy on the top (a) and 
bottom (b) surface layers, respectively. The red arrow 
contours denote the spin polarization vector of states at 
Ep, with the remaining energy contours representing an 
energy range of ±25 meV about the Fermi energy. We 
find that S r ~ —Sr and Sr is locked almost nor- 

-/C|| 

mal to the in-plane momentum fc||, which results from a 
combination of the strong SOC and the inversion sym- 



metry breaking at the surface, the well-known Rashba 
effect. [36] This in turn results in the rotation of the spin 
orientation around the Fermi surface (spin chirality) as 
reported for other known TI materials. 



IV. CONCLUSIONS 

In conclusion, we have investigated the evolution of the 
band structure and of the topological invariant in HgSe 
under application of [001] and [110] uniaxial strain. We 
predict that HgSe is a strong TI (i/q =1) in the com- 
pressive strain range, efo'oiKiiio]) < e[ooi]([iio]) < 0, where 
the (crystal field) band gap occurring away from F, dis- 
plays a non-monotonic strain behavior. The band gap 
vanishes at e[ooi]([iio])j where the system undergoes a TI 
o NI phase transition. We find that HgSe exhibits a 
similar behavior applying a [110] tensile strain. On the 
other hand, under [001] tensile strain HgSe remains a NI. 
Thus, these calculations demonstrate that HgSe can be 
tuned into a 3D topological insulator via proper strain 
engineering. 
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